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Abstract 

We give a simple algorithm to find a spanning tree that simultaneously approxi- 
mates a shortest-path tree and a minimum spanning tree. The algorithm provides a 
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Q , continuous trade-off: given the two trees and a 7 > 0, the algorithm returns a span- 

ning tree in which the distance between any vertex and the root of the shortest-path 
tree is at most 1 + y/2') times the shortest-path distance, and yet the total weight of 
the tree is at most 1 + \/2/7 times the weight of a minimum spanning tree, 
^-j. \ Our algorithm runs in linear time and obtains the best-possible trade-off. It can 

be implemented on a CREW PRAM to run in logarithmic time using one processor 
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per vertex. 
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1 Introduction 



A minimum spanning tree of an edge-weighted graph is a spanning tree of the graph of 
minimum total edge weight. A shortest-path tree rooted at a vertex r is a spanning tree 
such that for any vertex v, the distance between r and v is the same as in the graph. 
Minimum spanning trees and shortest-path trees are fundamental structures in the 



study of graph algorithms [11], [15], [19|]; fast algorithms for finding each are known 
]12 , |13| . Typically, the edge- weighted graph G represents a feasible network. Each vertex 
represents a site. The goal is to install links between pairs of sites so that signals can be 
routed in the resulting network. Each edge of G represents a link that can be installed. 
The cost of the edge reflects both the cost to install the link and the cost (e.g., time) for a 
signal to traverse the link once the link is installed. A minimum spanning tree represents 
the least costly set of links to install so that all sites are directly or indirectly connected, 
while a shortest-path tree represents the set of links to install so that for each site, the cost 
for a signal to be sent between the site and the root of the tree is as small as possible. 

The goal of a minimum spanning tree is minimum weight, whereas the goal of a shortest- 
path tree is to preserve distances from the root. We show that a single tree can approxi- 
mately achieve both goals. That is, the cost to install a set of links so that every site has 
a short path to the root is only slightly more than the cost just to connect all sites. 

In the graph in Figure J], the weight of the shortest-path tree is much more than the 
weight of a minimum spanning tree. Conversely, in the minimum spanning tree, the distance 
between the root and one of the vertices is much larger than the corresponding shortest- 
path distance. Nonetheless, there is a tree which nearly preserves distances from the root 
and yet weighs only a little more than the minimum spanning tree. We call such a tree a 
Light Approximate Shortest-path Tree (LAST). The main result of this paper is that such 
trees exist in all graphs and can be found efficiently. 

Let G = (V, E) be a graph with non-negative edge weights and a root vertex r. Let G 
have n vertices and m edges. Let u>(e) be the weight of edge e G E. The distance D G {u, v) 
between vertices u and v in G is the minimum weight of any path in G between them. 

Definition 1 For a > 1 and (3 > 1, a spanning tree T of G meeting the following two 
requirements is called an (a,/3)-LAST rooted at r. 

• (Distance) For every vertex v, the distance between r and v in T is at most a times 
the shortest distance from r to v in G. 

• (Weight) The weight ofT is at most (3 times the weight of a minimum spanning tree 
ofG. 



Theorem 1 (Section |3|) Let G be a graph with non-negative edge weights; letr be a vertex 
of G; let a > 1 and [3 > 1 + -^j. Then G contains an (a, f3)-LAST rooted at r. The LAST 
can be computed in linear time given a minimum spanning tree and a shortest-path tree, 
and in Oim + nlogn) time otherwise. 
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(a) Euclidean graph 
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(b) Min. spanning tree 
with distance blow-up 




(c) Heavy shortest-path tree 



(d) Light, approximate 
shortest-path tree 



Figure 1: Approximating both a minimum spanning tree and a shortest-path tree. 
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Note that there is a trade-off between the approximations of the two trees. The trade-off 
is the best possible: 

Theorem 2 (Section Fix a > 1 and 1 < (3 < 1 + -^-r. There exists a planar graph G 
with a vertex r such that G contains no (a, (3)-LAST rooted at r. Deciding whether a given 
graph contains an (a, (3) -LAST rooted at a given vertex is NP-complete. 

Note that for (3 = 1, the problem is to find a minimum spanning tree that best approx- 
imates a shortest-path tree. It follows from Theorem [| that this is NP-complete. When 
a = 1, the problem is to find a minimum- weight shortest-path tree. This problem can be 
solved in linear time, even in directed graphs: 

Theorem 3 (Section |5|) Given any shortest-path tree of a directed or undirected graph 
rooted at a given vertex, a minimum-weight shortest-path tree can be found in linear time. 

Finally, LAST's can also be found quickly in parallel, given a minimum spanning tree 
and shortest-path tree (or approximations thereof, see Section |3T^ ): 

Theorem 4 (Section |6|) Given a > 1, a minimum spanning tree, and a shortest-path 
tree, an (a, 1 + -^-r)-LAST can be found by n processors in O(logn) time on a CREW 
PRAM. 

2 Related Work 

Trees realizing tradeoffs between weight and distance requirements were first studied by 
Bharath-Kumar and Jaffe ||. The authors' weight requirement was the same as ours, but 
their distance requirement was that the sum of the distances from the root to each vertex 
should be at most (3 times the minimum possible sum. They showed the weaker tradeoff 
that the desired tree exists if a(3 > 0(n). 

Awerbuch, Baratz and Peleg ||, motivated by applications in broadcast-network design, 
made a fundamental contribution by showing that every graph has a shallow-light tree - 
a tree of diameter at most a constant times the diameter of G and of weight at most a 
constant times the weight of the minimum spanning tree. Our algorithm is modification 
of their algorithm. Cong, Kahng, Robins, Sarrafzadeh and Wong [^, |8|, |9|, motivated by 
applications in VLSI-circuit design, improve the constants in the construction of || and 
consider variations bounding the radius of the tree instead of the diameter. Recently and 
independently, Awerbuch, Baratz and Peleg 0, modified the algorithm from [[J. They 
obtained the same algorithm as in |J but a stronger analysis, proving that the algorithm 
computes an (a, 1 + ^-j-)-LAST. Their algorithm takes 0(m + n log n) time. Our algorithm 
achieves a strictly stronger distance requirement than the above algorithms. 

Considerable research has been done on finding spanners of small size and weight in 
arbitrary graphs 0, ||], |l7] and in Euclidean graphs (graphs induced by points in the plane) 
H [| [16|, ^(J . A t-spanner is a low- weight subgraph G' of G such that for any two vertices, 
the distance between them in G' is at most t times the distance in G. It is known that there 
are graphs that do not have constant-spanners of net weight bounded by a constant times 
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the weight of the minimum spanning tree. Awerbuch, Baratz and Peleg [|3| also consider 
light trees that have low average distance-blowup on all non-tree edges. References to most 
of the work on graph spanners may be found in the paper by Chandra, Das, Narasimhan 
and Soares ||. 

We can reduce the problem of finding an (a, 1 + ^-j-)-LAST to the problem of finding 
an a-spanner of weight at most + ^zf) times the minimum spanning-tree weight in a 
planar graph. An algorithm achieving the latter is given in This gives an alternate 
(but less efficient) method of for finding an (a, 1 + ^-j-)-LAST. 



3 The Algorithm 

The algorithm is given an a > 1, a minimum spanning tree, and a shortest-path tree rooted 
at a vertex r. It returns an (a, 1 H — ^M-LAST rooted at r. 

The basic idea of the algorithm is to traverse the minimum spanning tree, maintaining 
a current tree, and checking each vertex when it is encountered to ensure that the distance 
requirement for that vertex is met in the current tree. If it is not met, the edges of the 
shortest path between the vertex and the root are added into the current tree. Other edges 
are discarded so that a tree structure is maintained. 

After all vertices have been checked and paths added as necessary, the remaining tree 
is the desired LAST. The final tree is not too heavy because a shortest-path is only added 
if the path that it replaces is heavier by a factor of a > 1. This allows a charging argument 
bounding the net weight of the added paths. 



3.1 Relaxation. 

The tree is maintained by keeping a parent pointer p[v] for each non-root vertex v. To 
avoid recomputing shortest-path distances when a path is added, the algorithm maintains 
a distance estimate d[v] for each vertex v. This distance estimate, which is an upper bound 
on the true distance in the current tree, is used in deciding whether to add a path to the 
vertex. The parent pointers and distance estimates are initialized and maintained as in 



T(| Section 25.1]: 
INITIALIZE0 

Initialize distance estimates, parent pointers. 

1 for each non-root vertex v do p[v] <— nil; d[v] <— oo 

2 d[r] ^ 

Relax(m, v) 

Check for shorter path to v through (u, v). 

1 if d[v] > d[u] + w(u,v) 

2 then d[v] <— d[u] + w(u, v) 

3 p[v] <— u 

After executing Initialize, the algorithm builds and updates the tree and maintains the 
distance estimates by a sequence of calls to Relax. The important invariant maintained 
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by Relax is that the edges {(p[u],u) : d[v] ^ 00} form a tree, with d[v] an upper bound 
on the distance between the root and v in the tree. 

3.2 The algorithm as a sequence of relaxations. 

The algorithm performs a depth-first search of the minimum spanning tree starting at the 
root. For a tree, a depth-first search is simply an edge- by-edge walk from the root vertex 
through the vertices of the tree. Each edge is traversed twice: once in each direction. At 
any time in the search, the sequence of edges traversed so far forms a walk (a non-simple 
path) through the visited vertices. The walk starts at the root and ends at a vertex that 
we call the current vertex. We also say the algorithm is visiting this vertex. 

The relaxations done by the algorithm are of two kinds. The first kind adds shortest 
paths. The first time vertex v is visited, if d[v] exceeds a times the distance from the root 
to v in the shortest-path tree, then the edges of the shortest path are relaxed as needed to 
lower d[v] to the shortest-path distance. 

The second kind extends or modifies the current tree to use a minimum-spanning-tree 
edge if it is useful. Specifically, when an edge (u,v) is traversed from u to v, Relax(m, v) 
is called. This guarantees inductively that d[v] is bounded by the weight of the shortest 
path from the root r to vertex v' plus the weight of the minimum-spanning-tree path from 
v' to v, where v' is the most recent vertex to have its shortest path added. This invariant 
is what allows the weight of the added paths to be bounded. 

When the depth-first search finishes, the current tree is the desired LAST. The full 
algorithm is given in Figure 0. 

3.3 A sample execution. 

Figure ^| shows a sample execution of the algorithm with a = 2 on the graph given in 
Frame (a). Frames (b) and (c) give, respectively, a minimum spanning tree (of weight 60) 
and a shortest-path tree. 

Initially all parent pointers are nil and each d[v] is infinite. The depth-first search of the 
minimum spanning tree visits the vertices in increasing order of their labels and traverses 
the edges of the minimum-spanning tree in the following order: 

(1,2), (2, 3), (3, 4), (4, 5), (5, 4), (4, 6), (6, 7), (7, 6), (6, 4), (4, 3), (3, 2), (2, 1), (1, 8), (8, 1) 

Recall that when an edge is traversed, it is relaxed. When a vertex is visited, if its current 
distance estimate is not small enough to guarantee the distance requirement then the edges 
on the shortest path to the vertex are relaxed, modifying the current tree. 

Frame (d) shows the state of the algorithm just after vertex 5 has been visited: the 
edges (1,2), (2,3), (3,4), and (4,5) were relaxed as they were traversed. Because d[5] was 
equal to 40 (more than twice the shortest-path distance) when vertex v was visited, edge 
(1,5) was relaxed, changing vertex 5's parent to vertex 1, and changing d[5] to 15. 

Frame (e) shows the state after vertex 7 — the next vertex to have its shortest path 
added — has been visited. Note that when edge (5,4) was traversed, from 5 to 4, its 
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FlND-LAST(T M ,T 5 ,r,a) 

Input: Min. spanning tree Tm, shortest-path tree Tg, vertex r, a > 1. 
Output: an (a, 1 + ^j)-LAST rooted at r. 

1 INITIALIZE0 

2 DFS(r) 

3 return tree T = {(f,p[f]) | v E V — {r}} 
DFS(w) 

Traverse the subtree of Tm rooted at u, relaxing edges as they are tra- 
versed, and adding paths from T$ as needed. 

1 if d[u] > a D Ts (r, u) 

2 then Add-Path(w) 

3 for each child v of u in Tm 

4 do Relax(m,-u) 

5 DFS(v) 

6 Relax(w,-u) 

Add-Path(-u) 

Relax edges along path from r to v in T5. 

1 if d[v] > D Ts (r,v) 

2 then ADD-PATH(parent Ts (w)) 

3 Relax (parent Ts (v),v) 



Figure 2: Algorithm to compute a LAST. 

relaxation changed vertex 4's parent to vertex 5 and updated d[4] to reflect the new shorter 
path (1, 5), (5, 4). The algorithm then traversed and relaxed edges (4, 6) and (6, 7), bringing 
vertices 6 and 7 into the tree. When vertex 7 was encountered, its distance estimate (40) 
exceeded twice the shortest-path distance (15), so the edges on the shortest path (1, 8), (8, 7) 
to vertex 7 were relaxed in that order. This added these edges to the current tree and 
brought down the distance estimates of these vertices. 

Frame (f) shows the final state of the algorithm. The parent pointers give the final 
tree. Note that the relaxation of edge (7, 6) from 7 to 6, changed vertex 6's parent. This 
was the final change made to the tree. Subsequent relaxations made by the traversal had 
no effect. Remaining distance estimates were small enough to guarantee that the distance 
requirements were met, so that Add-Path was not called. 
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(a) Graph 



(b) Min. spanning tree 



(c) Shortest-path tree 




Figure 3: A sample execution of the algorithm. 
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3.4 Analysis of the algorithm. 

Next we prove that Find- Last (T^/, T$, r, a) returns an (a, 1 + ^j-)-LAST in linear time. 
Let T be the tree returned. 

Lemma 3.1 The distance between v and r in T is at most a times the shortest-path dis- 
tance. 

Proof. When a vertex v is visited, if d[v] exceeds a times the distance in the shortest-path 
tree then Add-Path is called, after which d[v] equals the shortest-path distance. In any 
case, after v is visited, d[v] is at most a times the shortest-path distance and subsequently 
never increases. On termination it bounds the distance in T. □ 

An amortized analysis establishes that the total weight of T is not too large. 

Lemma 3.2 The weight of T is at most (l + -^zi) times the minimum spanning-tree 
weight. 

Proof. Let v = r and let v±, v 2 , v k be the vertices that caused shortest paths to 
be added during the traversal, in the order they were encountered. When the shortest 
path from r to Vi (i > 1) was added, the net weight of the added edges was at most 
D Ts (r, v^. Also, the edges on the path to v^ consisting of the shortest path to v j_i followed 
by the path in the minimum spanning tree from to t>, had been relaxed in order, 
so that d[vi] < D Ts (r,Vi-i) + DT M (vi-i,Vi). The shortest path to Vi was added because 
aDT s (r, < d[vj\. Combining the inequalities, 

aD Ts (r,v,i) < D Ts (r,Vi-i) + Z> Tm (%_i, v t ). 

Summing over i bounds the net weight of the added paths: 

k k 

a ^Z D T s {r,Vi) < J2( D T 3 (r,v l - 1 ) + D Tu (vi-i,Vi)) 
i=i i=i 

and therefore 

k k 
(a-1) ^D Ts (r,Vi) <J2 D T M {vi-i,Vi). 

i=l i=l 

The DFS traversal traverses each edge exactly twice, and hence the sum on the right-hand 
side is at most twice the weight of Tm, i-e., 

k 

Y^D TM (vi„i,Vi) < 2w(T M ). 
i=i 

Hence the net weight of the added paths is less than -^-?w(Tm)- □ 
The following alternate proof of Lemma |3.2| may also be of interest . 

Alternate Proof of Lemma p. | . As the algorithm executes, define the potential function 
$ to be the distance estimate of the current vertex. When a shortest path of length p to 
the current vertex v is added, $ = d[v] > a p. Adding the path lowers d[v] to p, decreasing 
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$ by at least (a — 1) p. Hence the total weight of the added paths is bounded by the sum 
of the decrements to $ during the course of the algorithm, divided by a — 1. 

Since $ is initially and always non-negative, the sum of the decreases is at most the 
sum of the increments. $ increases only when the current vertex changes from some vertex 
u to a vertex v after the edge (u,v) was relaxed. This ensures that d[v] < d[u] + w(u,v) 
and that $ increases by at most w(u,v). Since each edge is traversed twice, the total of 
the increases to $ during the course of the algorithm is bounded by twice the weight of the 
minimum spanning tree. 

This establishes that the total weight of the added paths is bounded by times the 
weight of the minimum spanning tree. □ 

The running time is proportional to the number of relaxations. This is 0(n) because 
each edge in Tm or T$ is relaxed at most twice by DFS and at most once by Add- Path. If 
the shortest-path tree and the minimum spanning tree are not given, they can be computed 
in 0(m + nlogn) time [15, This establishes Theorem [TJ 



Observation 1: In metric graphs (complete graphs with edge weights satisfying the tri- 
angle inequality, such as Euclidean graphs) the shortest-path tree is trivial and can be 
found in 0(n) time. For Euclidean graphs (induced by points in the plane) the minimum 
spanning tree can be computed in O (nlogn) time |18| . In these cases the LAST can be 
found more quickly. 

Observation 2: If the algorithm is given an a- approximate shortest-path tree and a b- 
approximate minimum spanning tree, the tree returned by the algorithm will be an (aa, b + 
2b/ (a — 1))-LAST. If such trees can be found more quickly then a LAST can also be found 
more quickly. 

Observation 3: In the multiple-root variant, the distance requirement is that in the final 
tree (or forest) the distance between each vertex and its nearest root should be at most a 
times the distance to any root in the original graph. This variant can be easily reduced to 
the original problem by adding an artificial root at distance from the multiple roots. 



4 Optimality of the Algorithm 

Next we show that the algorithm is optimal in the following sense. Fix a > 1 and 1 < /3 < 
1 + -^j. There is a planar graph not containing an (a,/3)-LAST rooted at a particular 
vertex. Further, it is NP-complete to decide whether a given graph contains an (a, /3)-LAST 
from a given root. 

4.1 Non-existence of L AST's when 8 < 1 + -^r. 

" a— 1 

Lemma 4.1 If a > 1 and 1 < (3 < 1 + then there exists a planar graph containing no 
(a, (3)-LAST rooted at a particular vertex. 

Proof. The graph is shown in Figure |j. The structure of the graph is as follows. The root 
r is connected to a central vertex c by a path of weight A of edges of weight some small 
5. The central vertex is connected through similar paths of weight B to the i leaves. The 
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Figure 4: A graph with no (a, /3)-LAST for 8 < 1 + ^y, (A = a + 1, B = a + e- l and 
C = 2). 

root is connected to each leaf with an edge of weight C. Let A — a+l, B — a + e — 1 
and C = 2, where e is an arbitrarily small constant. For small enough 5, the minimum 
spanning tree is formed by using all edges except those of weight C. Notice that this graph 
is planar. 

Consider the paths from the root to any leaf. The shortest path is the direct edge of 
weight 2. Any other path weighs more than 2a: the path through the center vertex weighs 
A + B = 2a + e; any path through another leaf weighs at least 2 + 2B = 2(a + e). This 
means that in any (a,/5)-LAST all £ edges of weight 2 are present. In addition, all but £ 
of the remaining edges are present. Therefore the weight of any (a,/3)-LAST is at least 
2£ + Tm~ £S, where Tm = (a + l)+£(a — 1 + e) is the weight of the minimum spanning tree. 
Hence the ratio of the weight of the (a, 8)-LAST to the weight of the minimum spanning 
tree is at least 

n *(2 -J) 

a + 1 + i(a - 1 + e) 

If P < 1 + ~~[) then the above exceeds (3 for sufficiently small e and 5 and sufficiently large 

e.u 

4.2 NP-completeness of LAST queries. 

Next we show that for any fixed a > 1 and 1 < f3 < 1 + is NP-hard to decide whether 
a given graph contains an (a, /3)-LAST rooted at a given vertex. Thus, it is unlikely that 
a polynomial-time algorithm exists for finding (a, /3)-LAST's when /3 < 1 + 

Clearly the problem is in NP. The proof of NP-hardness is in two parts. We first 
show NP-hardness for /3 — 1 and fixed a > 1. We then reduce this problem to the fixed 
8 < 1 + -^t case. 

' a— 1 

Lemma 4.2 For /ixed a > 1, deciding the existence of an (a, 1)-LAST rooted at a given 
vertex of a given graph is NP-hard. 
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Figure 5: Reduction From 3-SAT. 
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Proof. The proof is by reduction from 3-SAT. Let F be a 3-SAT formula in conjunctive 
normal form — each clause consists of three literals from {x±, . . . ,x n } U {xi, . . . , x n }. We 
build a graph in which the (a, l)-LAST's correspond to satisfying assignments of F. 

A, B, D, E and W are constants to be determined later. The graph has a root vertex 
R, a vertex S, and a path connecting R to S of weight D consisting of edges small enough 
to ensure that the path is in any minimum spanning tree. 

For each pair of literals Xi and x~i, there are two vertices Xi and Xi, each having an edge 
to S of weight A. A path of weight E connects Xi and X{. This path is also constructed 
so as to be in any minimum spanning tree. 

For each clause Cj there is a vertex Cj with an edge to R of weight W. From Cj to each 
vertex corresponding to a literal in Cj there is an edge of weight B. 

This defines the graph. Observe that, provided < A < B < W, the minimum 
spanning trees are exactly characterized by the following. In any minimum spanning tree, 
the path from R to S and each path from Xi to X{ are present. For each variable Xj, 
exactly one of the two edges {(S, Xi), (S, Xi)} is present. For each clause Cj, exactly one 
edge of the form (JQ, Cj) or (Xi, Cj) for some i is present. No other edges are present. 

Next we use the distance requirement to ensure that any minimum spanning tree is an 
(a, 1)-LAST if and only if the edge to each clause vertex comes from some variable vertex 
Xi or Xi that has an edge in the minimum spanning tree directly to S. This is all that is 
needed, for then the (a, l)-LAST's will correspond to satisfying assignments in the original 
formula, and vice versa, as follows: for each variable Xi, choose the edge (S, Xi) iff Xi is 
true, otherwise choose the edge (S, Xi); for each clause Cj, choose the edge (Xi,Cj) (or 
(Xi,Cj)), where Xj (or Tj) is a variable (or negated variable) satisfying Cj. 

It suffices to choose A, B, D, E and W so that 

< A < B < W 

D + A + E < a min{A + D, B + W} 

D + A + B < a min{D + A + B,W} < D + A + E + B. 

To achieve this, let A = 1, B = a, D = 2a, E = (a - l)(2a + 1) and W = 1 + 2a + K □ 

Next we will reduce the (a, 1)-LAST problem to the (a, /3)-LAST problem, for any fixed 
a and j3 such that a > 1 and 1 < (3 < 1 + -^r. 

Proof. (Theorem |2|) Let G* be the graph for which we want to determine the existence 



of an (a, 1)-LAST rooted at a given vertex r*. By Lemma |4.1| , there exists a graph C 
with no (a, /3)-LAST rooted at some vertex r'. Assume without loss of generality that the 
minimum spanning tree of G* has weight 1 and the minimum spanning tree of C is of 
weight c (a constant to be determined later). Define the graph G to be the union of G* 
and G' by identifying r* and r' into a single root r. 

Let (3' be the minimum (3 such that G' has an (a,/3)-LAST. Define (3* analogously for 
G*. Take c = 

The weight of the minimum spanning tree in G is 1 + c; similarly the lightest tree in 
G meeting the distance requirement is of weight (3* + (3' c. Thus G has an (a, /3)-LAST 
iff /3* + (3' c < (3(1 + c). By our choice of c, this is equivalent to (3* < 1. Thus G has an 
(a,/?)-LAST iff G* has an (a, 1)-LAST. □ 
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5 Minimum- Weight Shortest-Path Trees 



Next we consider the case when a = 1, i.e., an (a, /3)-LAST is a shortest-path tree of weight 
at most (3 times the weight of the minimum spanning tree. In this case, no algorithm can 
guarantee any fixed /3 for all graphs. Instead, we show how to find a (l,/3)-LAST with 
minimum (3 in a given graph, i.e., a minimum- weight shortest-path tree. 

In fact, we solve a more general problem: finding a minimum- weight shortest-path tree 
in a rooted directed graph. The undirected case reduces to this case by the standard trick 
of replacing each undirected edge {u, v) by two new directed edges {u, v) and (v, u) of the 
same weight as the original edge. 

The directed problem reduces in turn to the problem of finding a minimum-weight 
branching in the shortest-path subgraph of the given directed graph. A branching is a 
directed spanning tree with all edges directed away from the root. The shortest-path 
subgraph is the spanning subgraph consisting of all directed edges (u, v) on some shortest 
path from the root, i.e., those for which Dc(r,u) + w{u,v) = Dc(r,v). It is easy to show 
that the shortest-path trees in a directed graph are exactly the branchings from the root in 
its shortest-path subgraph. Consequently, it suffices to find a minimum-weight branching 
in the shortest-path subgraph. 

A polynomial-time algorithm for finding a minimum-weight branching in any given 
graph is known ]13[ . However, a shortest-path subgraph of a non- negatively weighted graph 
has the property that any edge on a cycle has weight zero. This allows the following linear- 
time algorithm. First, identify the strongly-connected components in the subgraph induced 



by the edges of weight zero. This can be done in linear time |10[ . For each component not 
containing the root, choose the minimum- weight incoming edge and call the vertex with an 
incoming chosen edge the base vertex of the component. For the component containing the 
root, call the root vertex the base vertex. For each component, find a branching of weight 
zero edges rooted at the base in the subgraph induced by the component. Finally, return 
the chosen edges together with the edges of the components' branchings. 

This set of edges forms a branching: each non-root vertex has an incoming edge and 
there are no cycles. The branching is of minimum weight because in any branching every 
non-root component has at least one incoming edge. It is straightforward to implement the 
algorithm to run in 0(n + m) time. This proves Theorem |3[ 



6 Finding LAST's in Parallel 

Given a > 1, a minimum spanning tree, and a shortest-path tree, an (a, 1 H — ^-j-)-LAST 
can be found using n processors in 0(logn)-time. The model of computation we use is 
the Concurrent- Read, Exclusive-Write Parallel RAM, in which independent, synchronized 
parallel processors share a common memory WM. Multiple simultaneous accesses to the 



same memory location are allowed only if all of the accesses are read operations. 

The algorithm is as follows. Let C = (ei, e 2 , e 2n _2) be the (directed) edges of the 
walk through the graph implicit in the depth- first search of the minimum spanning tree, 
as in Section |3.2|. This tour can be constructed in O(logn) time by n processors using 
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standard techniques |L4|]. Let (iii,Ui + i) = ej. Using the terminology of Section j3T2| , after 



edge (ui,Ui + i) is traversed from m to Ui+i, vertex m+i is the current vertex. 

The parallel algorithm emulates the serial algorithm except that the distance estimates 
are looser in two ways. First, while a vertex may occur several times in C, the parallel 
algorithm treats each occurrence as a distinct vertex. Second, when a shortest path is 
added, only the distance estimate of the destination vertex is lowered. These looser distance 
estimates can be computed in parallel, but still suffice to imply the weight requirement. 

Let Dc(ui,Uj) denote Y?k=^ w ( e k), the distance from ui to Uj along C. Let m(i,j) be 
the relation 

i < j and DT s (r,Ui) + Dc(ui,Uj) > a DT 3 (r,Uj). 

The meaning of m(i,j) is the following. Suppose we modify the original algorithm to use 
the looser distance estimates. When the modified algorithm enountered a vertex Uj, if the 
algorithm had not added shortest paths to any of the vertices tij+i, Wj+2, Uj-i, then it 
would add the shortest path from the root to Uj. Thus, if the modified algorithm adds a 
shortest path to a vertex m, then the next shortest path it adds will be to vertex Uk, where 
k = min{j : m(i,j)}. 

The parallel algorithm will emulate the modified algorithm. Define J(i) = min{j : 
m(i,j)}. The parallel algorithm will compute the function J and then add shortest paths 
to vertices in the set 

S = {Ui, Mj(J(l)), Uj(J(....J(l). ..))}• 

Once J has been computed, S can be computed by n processors in O(logn) times on a 
CREW PRAM using standard techniques JL4J. Once S has been computed, the set S* of 



ancestors of S in the shortest-path tree can also be computed by n processors in O(logn) 



time using tree contraction techniques jlj] . The final tree is formed by each non-root vertex 
choosing as its parent either the parent in the shortest-path tree (if the vertex is in S*) or 
the parent in the minimum spanning tree (otherwise). It can easily be shown that every 
vertex has a path to the root using this set of n — 1 edges, so that they do indeed form a 
tree. 

It remains to compute J(i). First, note that m(i,j) is monotone in i for fixed j. 
Lemma 6.1 If i' < i and m(i,j) is true then m(i',j) is true. 

Proof. If m(i,j) is true then D Ts (r,Ui) + D c {ui,Uj) > a D Ts {r,Uj). For any i' < i, 
Dc{ui,Uj) < Doiu'^Uj). The shortest path from r to % is no longer than any other path 
from i to r in the graph, and hence DT s (r, ui) < Dt s (t, U{>) + Dciu'^ U{). Combining these 
inequalities, we get 

D Ts {r,u v ) + Dciu'^Uj) = D Ts {r,u v ) + D c {u' i ,u i ) + Dc{ui,Uj) 

> D Ts (r,Ui) + Dc(ui,Uj) 

> aD Ts (r,Uj) 

Hence m(i',j) is true by definition. □ 

The function J can be computed efficiently because of this monotonicity property: 
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Lemma 6.2 Suppose m(i,j) implies m(i f ,j) for 0<i'<i<n,0<j<n. Then the 
function J(i) = min{j : m(i,j)} can be computed in O(logn) time by n processors on a 
CREW PRAM. 



Proof. Define I (J) = max{i : m(i,j)}. For each j, compute using binary search. 
Define I* (J) = max.{I(j') : 1 < j' < j}. Compute function I* from function / using a 
standard prefix-maxima computation. Finally, define J'(i) = min{j : > i}- Compute 
function J', again using binary search, from monotone function J*. (See Figure ^|.) Each 
of these computations can be done by n processors in O(logn) time on a CREW PRAM 
using standard techniques [0 . 



J(4)=5 



m(i,j) 



Key: 



Q m(ij) true 
fj m(i,j) false 



- I(j) 
| J(i) 



1 1 1 1 1 1 — \+- 

1 2 3 4 5 6 7 

I(7)=3 



i \ i 



10 

I*(9)=6 



Figure 6: Computing J from m to find a LAST in parallel. 

We prove that J'(i) = J(i) for each i. The proof is in two steps. 

1. J'(i) = min{j : > i} = min{j : > i}, i.e., the smallest j such that 
exceeds i is equal to the smallest j such that exceeds %. This is because the latter 
depends only on the maxima of / — those j such that > I(j') for all j' < j. 

2. min{j : > i} = min{j : m(i,j)} because > i is equivalent to m(i,j) by the 
monotonicity property of m. 



The analyses of Lemmas |3.1| and |3.2| can easily be adapted to prove that the final tree 
produced by the parallel algorithm is an (a, 1 + ^-j-)-LAST. This establishes Theorem f| 
- an {a, 1 + ^)-LAST can be computed by n processors in O(logn) time on a CREW 
PRAM. 
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7 Conclusions 



Every graph contains trees that offer a continuous tradeoff between minimum spanning trees 
and shortest-path trees. Trees achieving the optimal trade-off can be found in (sequential) 
linear time or in logarithmic time by linearly many processors. 

In the following cases, is it possible to obtain a better trade-off? 



• In Euclidean graphs (i.e., points in the plane). Note that the proof of Lemma |3]2 
requires only that the algorithm walk around the graph from the root visiting every 
vertex once, i.e., that the algorithm traverse a Traveling Salesman path starting at 
the root. In Euclidean graphs, perhaps such a path of weight at most (2 — e) times the 
minimum spanning-tree weight always exists and can be found in polynomial time. 

• If the distance requirement is replaced by the requirement that the sum of distances 
from the root is within a times the minimum possible? 

• If the root is not fixed? This would correspond to the problem of installing a low-cost 
network and choosing a root site so that distances from the root are near-minimum. 

Clearly any (a,/3)-LAST also meets these looser requirements, but our lower bounds no 
longer show that the trade-off is optimal. 

For directed graphs, it is easy to show that for any fixed a and (3, (a, /3)-LAST's may 
not exist and that finding the minimum (3 such that an (a,/3)-LAST exists is NP-hard. 
Can one approximate this minimum (31 
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